Formation of Nuclear "Pasta" in Supernovae 
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In supernova cores, nuclear "pasta" phases such as triangular lattice of rod-like nuclei and layered 
structure of slab- like nuclei are considered to exist. However, it is still unclear whether or not they 
are actually formed in collapsing supernova cores. Using ab-initio numerical simulations called the 
Quantum Molecular Dynamics (QMD), we here solve this problem by demonstrating that a lattice 
of rod- like nuclei is formed from a bcc lattice by compression. We also find that, in the transition 
process, the system undergoes zigzag configuration of elongated nuclei, which are formed by a fusion 
of two original spherical nuclei. 

PACS numbers: 26.50.+X, 21.65.-f, 02.70.Ns, 97.60.Bw 



The mechanism of collapse driven supernova explosions 
[H has been a long-standing mystery Matter in super- 
nova cores is also yet to be understood and these issues 
are closely connected. In the initial stage of the explo- 
sions, the collapsing iron core experiences an adiabatic 
compression, which leads to an increase of the central 
density from ~ 10^ g cm~^ to around the normal nu- 
clear density ~ 3 x 10^^ g cm~^ (corresponding to the 
nucleon number density po = 0.165 fm~^) just before 
the star rebounds. Thereby, the Coulomb repulsion be- 
tween protons in nuclei, which tends to make a nucleus 
deform, becomes comparable to the surface tension of 
the nuclei, which favors a spherical nucleus. The pasta 
phases (phases consisting of "spaghetti" -like columnar 
nuclei and of "lasagna"-like planar nuclei, etc.) [1, 0, HI 
are thus expected to be formed in the inner cores dur- 
ing the collapse of stars and amount to more than 20% 
of the total mass of the cores just before the bounce [6]. 
Since the coherent scattering of neutrinos is very differ- 
ent between pasta nuclei and spherical ones 0, S] , pasta 
phases can have a significant influence on the dynamics 
of explosions 0, Q . Thus there is a growing interest in 
pasta phases [9| . 

However, the above speculation that the pasta phases 
are formed in collapsing cores is based on p hase diagrams 
of the equilibrium state (e.g., Refs. [T^,!!^,!!!!,!!^, 14] for 
non-zero temperatures) or static and perturbative calcu- 
lations [iBI, ll6| . Since formation of the pasta phases from 
a bcc lattice of spherical nuclei is accompanied by dy- 
namical and drastic changes of the nuclear structure, the 
fundamental question whether or not the pasta phases 
are formed in supernova cores is still open, and an ab- 
initio approach is called for. In the present work, we 
answer this question using called the quantum molecular 
dynamics (QMD) 171, ll8|. QMD can properly incorpo- 



rate the thermal fluctuations and is a powerful and suit- 
able approach [19] to describe the deformation of nuclei 
in the present problem as has been exploited for studying 
the equilibrium phase diagram 12|, 2Q|] and the transition 
dynamics between the pasta phases [21| . 

We use the QMD Hamiltonian of Ref. with the 
standard medium-equation-of-state parameter set [22|]. 
In our simulations, we consider a system with protons, 
neutrons, and charge-neutralizing electrons in a cubic 
box with periodic boundary condition. The electrons 
are relativistic and degenerate, and can be well approx- 
imated as a uniform background 0, ^ . We calculate 
the Coulomb interaction by the Ewald sum. Our simula- 
tions are carried out for the proton fraction x c^i 0.39 and 
0.49 using different initial conditions of the bcc lattice. 
In the following, we focus on the results for x 0.39 (the 
qualitative results are the same for x 0.49). 

The initial condition oi x 0.39 is obtained in the 
following way. We first prepare an isolated nucleus of 
^^^Pb at zero temperature and make a unit cell of the 
bcc lattice by setting two copies of the nuclei in a box 
with the periodic boundary condition. Using the Nose- 
Hoover thermostat for momentum-dependent potentials 
[l2| . we then equilibrate this sample at the temperature 
T = 1 MeV. We combine eight replica of this bcc unit 
cell to make a sample of the total number of nucleons 
N = 3328 (with 1312 protons and 2016 neutrons) at 
the nucleon number density p = 0.1 5po (the box size 
L = 51.23 fm) and equilibrate at T = 1 MeV for 7800 
fm/c using the Nose- Hoover thermostat and further relax 
for 4100 fm/c without the thermostat. Equilibrating 
this sample at different temperatures, we also prepare 
initial conditions for various temperatures. 

We simulate the compression of the bcc phase of spher- 
ical nuclei in the collapse. Starting from the above ini- 
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tial conditions, we increase the density by changing the 
box size L (the particle positions are rescaled at the 
same time). Here the average rate of the compression 
is < O{10~^) po/{fm/c) yielding the time scale of ^ 10^ 
fm/c to reach the typical density region of the phase with 
rod-like nuclei. This is much larger than the time scale 
of the change of nuclear shape (e.g., ~ 1000 fm/c for 
nuclear fission) and thus the dynamics observed in our 
simulation would be determined by the intrinsic physical 
properties of the system, not by the density change ap- 
plied externally. We perform adiabatic compression and 
isothermal compression at various temperatures [2^. In 
all the cases, we observe the formation of rod-like nuclei; 
here we show some typical examples in which we obtain 
a clear lattice structure of the rod-like nuclei. 

Figure [T] shows the snapshots of the formation process 
of the pasta phase in adiabatic compression. Here, we 
start from the initial condition of x = 0.39 and T = 0.25 
MeV {t = fm/c). At t 0^ 57080 fm/c and p 0^ 0.243po 
[Fig. [IJc)], the first pair of two nearest-neighbor nuclei 
start to touch and fuse (dotted circle), and then form 
an elongated nucleus [see, e.g., Fig.IlJd)]. After multiple 
pairs of nuclei become such elongated nuclei, we observe 
zigzag structure as shown in Fig. (H^d). These elongated 
nuclei start to stick together at t 59000 fm/c and p = 
0.246po, and all the nuclei fuse to form rod-like nuclei at 
t < 72700 fm/c and p < 0.267po. At t = 76570 fm/c and 
p = 0.275po7 we stop the compression (T :^ 0.5 MeV in 
this stage). We then relax the system microcanonically 
for 65510 fm/c and further relax at T = 1 MeV for :^ 
10000 fm/c (after that we take back the temperature to 
:^ 0.5 MeV and relax for 20000 fm/c; for the increase 
and decrease of the temperature between 0.5 and 1 MeV, 
we take :^ 5000 and 15000 fm/c, respectively). Finally, 
we obtain a triangular lattice of rod-like nuclei [Figs.[Hh- 
1) and (h-2)]. Remarkable point is that, in the middle of 
the transition process, pair of spherical nuclei get closer 
to fuse in a way such that the resulting elongated nuclei 
take a zigzag configuration (hereafter, we call it "zigzag 
pairing" ) and they further connect to form wavy rod-like 
nuclei. This feature is observed in all the other cases 
in which we obtain a clear lattice structure of rod-like 
nuclei (among them, the highest temperature case is the 
adiabatic compression from T = 1.5 MeV), and the above 
scenario of the transition process is qualitatively the same 
also for those cases. It is very different from a generally 
accepted picture (see, e.g., p. 462 of Ref. [15]) that all 
the nuclei elongate in the same direction along the global 
axis of the resulting rod-like nuclei and they join up to 
form straight rod-like nuclei. 

It has been regarded so far that the formation of the 
pasta phases in supernova cores is triggered by the fission 
instability with respect to the quadrupolar deformation 
of spherical nuclei [lH. To examine this point, we in- 
vestigate the variance of the radius of each nuclei over 
the solid angle and the area-averaged mean curvature of 




FIG. 1: (Color online) Snapshots (two-dimensional projec- 
tion) of the transition process from the bcc lattice of spher- 
ical nuclei to the pasta phase with rod-like nuclei. The red 
particles show protons and the green ones neutrons. In panels 
(a)-(g) and (h-1), nucleons in a limited region relevant for the 
two rod-like nuclei in the final state (h) [surrounded by the 
dotted lines in panel (h-2)] corresponding to eight nuclei in 
the initial condition (a) are shown for visibility. The vertices 
of the dashed lines in panels (a) and (d) show the equilibrium 
positions of nuclei in the bcc lattice and their positions in the 
direction of the line of sight are indicated by the size of the 
circles: vertices with a large circle, with a small circle, and 
those without a circle are in the first, second, and third lat- 
tice plane, respectively. The dotted circle in panel (c) show 
the first pair of nuclei start to touch. The solid lines in panel 
(d) represent the direction of the two elongated nuclei: they 
take zigzag configuration. In the final state [(h-1) and (h-2)], 
almost perfect triangular lattice of rod-like nuclei is obtained. 
The box sizes are rescaled to be equal in the figures. 

the nuclear surface; neither of these quantities show a 
significant increase throughout the compression process 
until the nuclei start to touch [see also Fig. [T][b)]. This 
means that, before nuclei deform to be elongated due to 
the fission instability, they stick together keeping their 
spherical shape. This is consistent with the result of Ref. 
[28] (see also Ref. [29]), which shows that, within an in- 
compressible liquid-drop model, the condition for the fis- 
sion instability is not satisfied if one takes account of the 
background electrons. 

To understand the mechanism which triggers the for- 
mation of the pasta phases, we calculate the distance Tnn 
between centers of mass (for protons) of nearest-neighbor 
nuclei in the compression process of our simulations |23| • 
Figure [2] shows the mean value (rnn) and the standard 
deviation Arnn of rnn calculated for 16 nuclei in the sim- 
ulation box. Let us first focus on the result obtained for 
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FIG. 2: (Color online) Distance rnn between centers of mass of 
nearest- neighbor nuclei in the compression process. The blue 
Une shows the mean value (rnn) and the blue filled region 
shows the standard deviation Arnn of rnn normalized by 
(a value of rnn for a perfect bcc lattice without displacement of 
nuclei). Panel (a) is for the isothermal compression at T = 
MeV and panel (b) is for the adiabatic compression starting 
from T — 0.25 MeV. The first pair of nuclei touch at a density 
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in the red filled region (corresponding to r 
both cases) whose width shows its uncertainty. 
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the isothermal compression at T = MeV [Fig. [2](a)], 
where we can see a clear signature due to the absence of 
thermal fluctuations. We note that there is a significant 
increase of Arnn just before the nuclei start to connect. 
This shows that, before pairs of nuclei touch, nuclei dis- 
place from the equilibrium position of the bcc lattice and 
the bcc structure is spontaneously broken. As a result, 
each pair of nuclei approach to fuse and then form an 
elongated nuclei arranged in a zigzag configuration. At 
non-zero temperatures, thermal fluctuation smears the 
above signature and it would also assist triggering the 
formation of the pasta phases; however even in the case 
of the adiabatic compression starting at T = 0.25 MeV, 
we can still see a significant increase of Arnn just before 
the nuclei touch and fuse [Fig. [2l^b)]. 

Using a simplified model, we now examine our results 
of the QMD simulations. When nearest neighbor nu- 
clei are so close that the tails of their density profile 
overlap with each other, net attractive interaction be- 
tween these nuclei starts to act due to the interaction 
between nucleons in different nuclei in the overlapping 
surface region. We consider that this attraction between 
nuclei leads to the spontaneous breaking of the bcc struc- 
ture. In order to examine this hypothesis, we consider 
a minimal model in which each nucleus is treated as a 
point charged particle interacting through the Coulomb 
potential and the potential of the Woods- Saxon form: 
V{r) = Vb {1 + exp [(r — i?)/a]}~^, which describes the 



•/ 



• 



(c) L=86.91 fm, t=252650 fm/c (d) L=86.91 fm, t=347160 fm/c 




FIG. 3: (Color online) Snapshots of the transition process in 
the simulation of isothermal compression at T < 0.01 MeV 
using the simplified model. Starting from the bcc lattice of 
nuclei (a), the first pairing (b), zigzag pairing around the first 
pair (c), and finally a zigzag structure (d) are observed. The 
red particles show the centers of mass of nuclei and the nuclei 
within the distance less than 0.89ri*5 are connected by a blue 
line. Nuclei and connections within only two lattice planes 
normal to the line of sight are shown. 



finite size of nuclei and models the internuclear attraction 
when the nearest neighbor nuclei start to touch. Since 
nuclei start to connect before they are deformed, it is 
reasonable to treat a nucleus as a sphere and incorporate 
only its center-of-mass degree of freedom. Parameters R 
and a represent the interaction radius of a nucleus and 
range of the effective nuclear forces, respectively and we 
take Vo = -11.5 MeV, R = 17 fm, and a = 0.5 fm (3^. 
With this model, we carry out the isothermal (keeping 
the system cool at T < 0.01 MeV by the frictional relax- 
ation method throughout the simulation; this small but 
non-zero temperature is sufficient for symmetry break- 
ing) compression of the bcc lattice of 128 nuclei of ^^^Pb, 
which corresponds to 8 times larger system than that 
of the QMD simulations. In Fig. [3] we show the snap- 
shots of this simulation. Here we compress from 0.221po 
(L = 90.00 fm) at a rate of p < (9(10-^) po/{fm/c). At 
0.246po {L = 86.91 fm) the first pair of nuclei starts to get 
closer [Fig. [SJb)] and then we stop the compression and 
relax the system. We observe zigzag pairing around the 
first pair [Fig. El^c)] and finally we obtain a zigzag struc- 
ture [Fig. [31(d)]. This supports that the zigzag structures 
observed in the QMD simulations are not caused by a 
finite size effect of the simulation box. 

In conclusion, we have shown that an ordered struc- 
ture of rod-like nuclei can be formed by compressing a 
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bcc lattice of spherical nuclei. Our result establishes that 
the pasta phases can be formed in collapsing supernova 
cores. Unlike a generally accepted conjecture, we have 
observed that nuclei start to connect before they deform 
due to the fission instability. This spontaneous breaking 
of the bcc structure is due to an attraction between nuclei 
caused by the overlap of the tails of nucleon distribution 
of neighboring nuclei. We have also discovered that, in 
the transition process, the system takes a zigzag configu- 
ration of elongated nuclei, which are formed by a fusion 
of original two spherical nuclei. Since a drastic change of 
the neutrino transport in the pasta phases pointed out 
in Ref. [7] has been already shown [8] (see also Ref. [g]), 
it is very interesting to perform core collapse simulations 
incorporating this effect as the next step. 
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